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Optical differentiation and optical Hilbert transformation play important roles in communications, 
computing, information processing and signal analysis in optical domain offering huge bandwidth. 
Meanwhile, silicon-based photonic integrated circuit is one of the most promising candidates for all- optical 
signal processing due to its intrinsic advantages of low power consumption, compact footprint, ultra-high 
speed and compatibility with electronic integrated circuits. In this study, we analyze the interrelation 
between first-order optical differentiation and optical Hilbert transformation and then experimentally 
demonstrate a feasible integrated scheme which can simultaneously function as first-order optical 
differentiation and optical Hilbert transformation based on a single microdisk resonator. This finding may 
motivate the development of integrated optical signal processors. 

Optical signal processing has received significant attentions in the past years in an attempt to overcome the 
bandwidth and speed bottlenecks incurred in conventional electronics 1 . Many optical signal processing 
systems have been proposed and demonstrated over the years, including all- optical logics 2 " 5 , optical 
differentiation (OD) 6 " 11 and integration 12 " 16 , and optical Hilbert transformation (OHT) 1718 . Particularly, OD has 
been widely used in optical computing, optical communications, optical metrology, optical digital processing and 
analog processing, and optical sensing 6 " 11 , while OHT has been used in the applications where signal phase is of 
critical importance, such as optical phase shifter 1718 , optical single sideband modulation (SSBM) 19,20 , vector 
modulation 21 and frequency measurement 22 . It is very interesting to discover the relationship between OD and 
OHT since these two functions have some similarities. A theoretical analysis of the interrelation of OD and OHT 
was discussed in Ref.18, but no further experimental verification was carried out yet. To date, most of the 
functionalities of OD and OHT, were implemented on the all-fiber platform (such as long periodical grating 7 ) 
or silicon photonics platform (such as microring resonator 11 ). The all-fiber platform can offer huge operation 
bandwidth (such as sub-picosecond pulse 23 ), very low propagation loss, and near null of coupling loss. The silicon 
photonics platform can offer compact footprint, well integration capability, and compatibility with 
Complementary Metal Oxide Semiconductor (CMOS) technology. Both these two platforms show their own 
advantages in all optical signal processing. In terms of large-scale integration, silicon-based photonic integrated 
circuit shows better performances. Therefore, in this paper we attempt to experimentally demonstrate OD and 
OHT simultaneously based on an on-chip microdisk resonator (MDR). The MDR has a very compact footprint of 
-20 urn X 20 urn. 

The transfer function of first-order OD can be expressed as 1 

H 0 d((o) =j(co-co c ) 
while for an ideal OHT, the transfer function is defined by 17 " 19 
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where j = V— 1, and cl> c are the optical angular frequency and 
central angular frequency, respectively. From Eqs. (1) and (2), we 
can see that the phase responses of the first- order OD and the ideal 
OHT are the same, with a phase shift of n at the central frequency. 
However, the amplitude response of the first- order OD is a notch 
response (linearly varied with the frequency) whereas that of the 
ideal OHT is an all-pass response, which could not be realized in 
practice. In fact, when we attempt to obtain an optical filter utilized 
for OHT, the optical filter normally turns to be a notch filter. Assume 
that the amplitude response is linearly varied with the frequency near 
the notch frequency. Then the transfer function is first- order OD at 
low frequency operation. However, for high frequency operation, the 
amplitude response can be regarded as a constant. Thus the transfer 
function satisfies OHT. Therefore, for a given notch filter with a n 
phase shift at the notch, it is determined by the bandwidth of input 
signals whether implementing OD or OHT. 

To quantitively investigate the interrelation between the first- 
order OD and OHT, we introduce a deviation analysis during simu- 
lation. Assume that the input Gaussian pulse signal has a fixed full 
width at half maximum (FWHM) of 24.52 ps and a repetition fre- 
quency of 6 GHz, as shown in Fig. 1(a). We define a notch filter 
whose 3-dB bandwidth is tunable, as shown in Fig. 1(b) (blue solid 
line) and the notch filter has a phase shift from n/2 to — n/2 at the 
central frequency (shown in Fig. 1(b), red dot line). Make the input 
Gaussian pulse propagate through the notch filter with its central 
frequency aligned to the notch. Then the output waveform of the 
notch filter is calculated, which will be compared with ideal output 
waveforms of OD and OHT results, respectively. 



To compare the waveform accuracy quantitively, an average devi- 
ation parameter is defined as the mean absolute deviation of the 
output waveform power from the ideal OD or OHT waveform at 
certain duration, which is —166.67 ps in our simulation. Then we 
vary the 3-dB bandwidth of the notch filter and scan the deviation 
parameters. Figure 1(c) shows the deviation as a function of the 3-dB 
bandwidth of notch filter, where blue circle represents the deviation 
between the output waveform and ideal OD and red square repre- 
sents the deviation between the output waveform and ideal OHT, 
and the insets cl-c3 show the output waveforms (red solid line) of 
the notch filter, calculated ideal OD (green dot line) and calculated 
ideal OHT (blue dash line) with different 3-dB bandwidths of the 
notch filter. 

Moreover, we also investigate a variety of different input pulse 
shapes in the same way. Figures 2(a) -(c) respectively show the devi- 
ation as a function of the 3-dB bandwidth of notch filter when the 
corresponding input pulse shapes are Sech 2 , super-Gaussian and 
Hermite-Gaussian, where blue circles represent the deviation 
between the output waveform and ideal OD and red squares repres- 
ent the deviation between the output waveform and ideal OHT, and 
the insets al-cl show the output waveforms (red solid line) of the 
notch filter, calculated ideal OD (green dot line) and calculated ideal 
OHT (blue dash line) with some certain 3-dB bandwidths of the 
notch filter. 

One can see that, for a larger bandwidth of the notch filter, such as 
35 GHz, the output waveform accords well with the calculated ideal 
OD result, otherwise, for a small bandwidth, such as 5 GHz, it closely 
approximates the calculated ideal OHT result. In other words, for a 
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Figure 1 | Simulation results of the first-order OD and the OHT. (a) input Gaussian pulse with fixed FWHM of 24.52 ps; (b) the spectrum of the input 
Gaussian pulse (green dash line), the amplitude response of the defined filter (blue solid line) with 3-dB bandwidth tunable and the phase response of the 
defined filter (red dot line); (c) the deviations between the output waveform and the OD (blue circle) and the OHT (red square) as a function of the 
defined filter's 3-dB bandwidth, insets (cl-c3) show the output waveforms (red solid line) of the defined filter, ideal OD waveform (green dot line) 
and ideal OHT waveform (blue dash line) with different 3-dB bandwidths of the defined filter. 
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Figure 2 | Simulation results of other input pulse shapes, (a), (b) and (c) the deviations between the output waveform and the OD (blue circle) and the 
OHT (red square) as a function of the defined filter's 3-dB bandwidth with input of Sech 2 pulse, super-Gaussian pulse and Hermite-Gaussian pulse, 
respectively; insets (al-cl) show the corresponding output waveforms (red solid line) of the defined filter, ideal OD waveform (green dot line) and ideal 
OHT waveform (blue dash line) with certain 3-dB bandwidths of the defined filter. 



bandwidth-fixed notch filter with a phase shift of n at the central 
frequency, a wideband input signal (small pulsewidth) will result in 
an OHT result, and a narrowband input signal (large pulsewidth) will 
result in an OD result on the contrary. 

Results 

Device structure. As illustrated in Fig. 3(a), an MDR structure fabri- 
cated on commercial silicon-on-insulator (SOI) wafer consists of a 
disk and a bus waveguide. The top silicon thickness of the SOI wafer 
we used is 340 nm. To form silicon ridge waveguide, the upper 
silicon layer is etched downward 240 nm, thus the lower silicon 
layer is 100 nm. Figure 3(b) shows the scanning electron micros- 
cope (SEM) image of the fabricated silicon MDR. The radius of the 
disk is 10 |im, and the waveguide width and thickness of the bus 
waveguide are 500 nm and 240 nm, respectively. The gap between 



the disk and the bus waveguide is 150 nm. We use vertical grating 
coupling method to couple the fiber with the silicon MDR, and the 
grating coupler we designed has a period of 630 nm with duty cycle 
of 56%. The coupling loss of the grating coupler is measured to be 
5 dB for a single side. Figs. 3(c) and 3(d) show the SEM images of the 
grating coupler and zoom-in coupling region, respectively. 

The measured transmission spectrum of the fabricated MDR is 
illustrated in Fig. 4(a), and Fig. 4(b) shows the zoom-in notch region 
selected in our experiment, which has a low loss of 12 dB, a relatively 
high extinction ratio (ER) of 15 dB and a relatively high Q factor of 
— 13688. The resonant wavelength is around 1554.17 nm. All the 
spectra are measured by an optical spectrum analyzer (OSA, 
AQ6370B) with a resolution of 0.02 nm. Figure 4(c) shows the mea- 
sured phase response of the MDR. One can see a 7r-phase shift is 
obtained. The method of phase measurement can be found in Ref. 24. 
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Figure 3 | MDR design, (a) schematic figure of the MDR composed of a 
disk and a bus waveguide (vertical grating couplers not shown); (b) SEM 
image of the fabricated MDR, the scale bar represents 10 (am; (c) SEM 
image of the vertical grating coupler, the scale bar represents 30 um; (d) 
SEM image of the zoom-in coupling region between the disk and the 
straight waveguide, the scale bar represents 2 (am. 
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Figure 4 | Characteristics of the MDR. (a) measured spectrum of the 
fabricated MDR; (b) the zoom-in spectrum of the notch we employed; (c) 
measured phase response of the MDR. 
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Figure 5 | Schematic of the experimental setup. First, a continuous wave beam is emitted by a tunable laser source and is then modulated into Gaussian- 
like pulses by the tunable pulse generator module. Subsequently, the signal is coupled into the MDR by the designed grating coupler. Finally, the 
MDR output is amplified and observed. TLS: tunable laser source. PC: polarization controller. RF: radio frequency. MZM: Mach-Zehnder modulator. 
EDFA: erbium doped fiber amplifier. OTDL: optical tunable delay line. PM: phase modulator. SMF: single mode fiber. OSO: optical sampling 
oscilloscope. 



Experiment overview. The experimental setup is shown in Fig. 5. 
The configuration mainly consists of two parts, namely, the tunable 
pulse generator and the MDR. The configuration of tunable pulse 
generator can be used to generate a Gaussian-like pulse with both 
repetition frequency and pulsewidth tunable. The output waveform 
of the MDR is amplified by an erbium doped fiber amplifier (EDFA) 
to compensate the loss of the MDR and then analyzed through a 
communication optical sampling oscilloscope (OSO). 

Input Gaussian pulse. The tunable pulse generator is composed of a 
tunable laser source (TLS), two polarization controllers (PCs), a 
Mach-Zehnder modulator (MZM) to produce a direct current 
(DC) free pulse train, an EDFA to compensate the loss induced by 
the MZM, an optical tunable delay line (OTDL) used to synchronize 
the radio frequency (RF) applied on the MZM and phase modulator 
(PM), a PM to induce large chirp, a segment of single mode fiber 
(SMF) for chirp compression to generate a short pulse train 25 . By 
changing the input RF power along with the length of the SMF, pulses 
with different FWHMs can be obtained. (See Methods for detailed 
explanation). 

Experimental result. The measured input and the output temporal 
waveforms are depicted in Figs. 6(a)-(h), and the calculated 
waveforms of ideal OD and OHT are also shown for comparison. 
Figures 6(a), (c), (e) and (g) show the input temporal pulses with 
FWHMs of 19.57 ps, 24.06 ps, 31.86 ps and 42.53 ps, respectively, 
fitting by Gaussian pulses. And Figs. 6(b), (d), (f) and (h) are the 
corresponding measured output waveforms at the MDR output port. 
One can see that when the pulsewidth of Gaussian pulse is as narrow 
as 19.57 ps, the output waveform agrees well with the ideal OHT 



waveform rather than the ideal OD waveform. On the contrary, when 
the input pulsewidith increases to 42.53 ps, the output waveform 
accords well with the ideal OD waveform. To analyze the output 
waveforms, we also calculate the average deviation of the measured 
output waveform from the calculated theoretical OHT and OD in 
one certain pulse period. Figure 7 shows the calculated average 
deviations of input pulse to standard Gaussian pulse (green 
triangle), output waveform to OHT waveform (red square) and 
output waveform to OD waveform (blue circle), respectively. With 
the increase of input pulsewidth, the deviation of output from the 
calculated theoretical OHT waveform increases, while that of output 
from the calculated theoretical OD waveform decreases. This trend 
has a good agreement with the theoretical analysis in Fig. 1(c). 

Additionally, time-bandwidth product (TBP) is the product of the 
FWHM of a signal and its 3-dB spectral width. We calculated the 
TBP of the output waveforms (red square), ideal OHT waveform 
(blue triangle) and ideal OD waveforms (green circle) with different 
input pulsewidths, respectively, as shown in Fig. 8. One can infer that 
the output TBP is closer to that of ideal OHT with a smaller input 
pulsewidth, whereas, it approximates the TBP of ideal OD for a larger 
input pulsewidth. 

Similarly, we also investigate the average deviation using the mea- 
sured frequency response of fabricated MDR for different input pulse 
shapes. Figures 9(a)-(c) show the calculated average deviations from 
ideal OHT waveform (red square) and from ideal OD waveform 
(blue circle) as a function of input pulsewidth when injecting 
Sech 2 , super- Gaussian and Hermite-Gaussian waveforms respect- 
ively, and the corresponding insets show the calculated output wave- 
forms (red solid line) of the MDR, ideal OHT (blue dash line) and 
ideal OD (green dot line) with different input pulsewidths. One can 
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Figure 6 | Experimental results, (a), (c), (e) and (g) are the inputs; red solid lines: input pulses, blue dash lines: fitted Gaussian pulses, (b), (d), (f) and (h) 
are the outputs; red solid lines: output waveforms, blue dash lines: calculated OHT, green dot lines: calculated OD. 
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Figure 7 | Analysis of the temporal waveforms of the MDR. Green 
triangle: average deviation between the input pulse and the fitted Gaussian 
pulse, red square: average deviation between the output waveform and the 
calculated theoretical OHT waveform, blue circle: average deviation 
between the output waveform and the calculated theoretical OD 
waveform. 

see that the deviation from OHT becomes much larger than that 
from OD with the increase of input pulsewidth, no matter what kind 
of waveform shape is injected. This trend also has a good agreement 
with the theoretical analysis in Fig. 2. 

Discussion 

In this work, we observed both the OD and OHT in a single SOI 
MDR chip. First, we numerically analyzed the interrelation between 
OD and OHT, as illustrated in Figs. 1 and 2. And we found that, for a 
properly designed notch filter, it can simultaneously function as OD 
and OHT under different input conditions. Because of a relatively 
high Q factor and ~tt phase shift at the resonated frequency, a 
properly designed MDR was chosen to demonstrate the OD and 
OHT, and the measured spectrum of the fabricated MDR was shown 
in Fig. 4. Then we performed the experiment. As illustrated in Figs. 6 
and 7, we confirmed in experiments that a single SOI MDR can 
function as an OD or an OHT as the input pulsewidth varied. 
Because of limited bandwidth of the MDR, the output waveform 
deviates from the ideal OD result under a short pulse injection. 
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Figure 8 | TBP Analysis of the temporal waveforms of the MDR. Red 

square: calculated TBP of the output waveforms, blue triangle: TBP of the 
ideal OHT waveforms, green circle: TBP of the ideal OD waveforms. 

Additionally, due to practically gradual transition of amplitude res- 
ponse at the notch, the output waveform deviates from the ideal OHT 
result under a wide pulse injection. All the experimental results 
accord well with the simulated results. 

In the part of Gaussian pulse generation, some chirps are induced 
by the PM firstly and then are compensated by the SMF. The input 
Gaussian pulse may still carry small chirps due to the incomplete 
compensation. Thus we need to analyze the impact of chirped pulse 
during the OD or OHT. From theoretical analysis, we find that 
whether the input Gaussian pulse is chirped or not, the calculated 
OD waveform remains the same. For the OHT, the output waveform 
may have very tiny distortion with asymmetry if the chirp parameter 
of the input Gaussian pulse is small. In the experiment, once we 
change the input RF power, we then try to generate a shortest pulse 
with constantly optimizing the length of SMF. Thus, it is reasonable 
to deduce that the chirp of the input Gaussian pulse is very small, 
without severe impacts on the output results. 

In summary, we have demonstrated an integrated scheme which 
can simultaneously function as OD and OHT based on a single SOI 
MDR. Also we have experimentally verified the interrelation between 
the OD and OHT. The scheme shows the advantages of compactness, 
flexibility, versatility and expandability due to the advantage of sil- 
icon-based technology. 
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Figure 9 | Analysis of different input pulse shapes, (a), (b) and (c) respectively show the calculated average deviations of calculated output waveform to 
ideal OHT waveform (red square) and output waveform to ideal OD waveform (blue circle) as a function of input pulsewidth when injecting Sech 2 , 
super-Gaussian and Hermite-Gaussian pulse. The insets (al-a2), (bl-b2) and (cl-c2) are the corresponding calculated output waveforms (red solid line) 
of the MDR, ideal OD (green dot line) and ideal OHT (blue dash line) with different input pulsewidths. 
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Methods 

Input Gaussian pulse. We used a tunable laser source with a tuning resolution of 
0.01 nm, which enables us to precisely align the signal wavelength to the resonance 
peak of MDR that we employ. The continuous wave (CW) light from the TLS is first 
modulated by the MZM with sinusoidal RF signal driven. To obtain a DC free pulse 
from the MZM, the driving parameters of the MZM should be carefully adjusted. 
Utilizing the nonlinear transmission curve of the MZM, a return-to-zero (RZ) pulse 
will be obtained when the MZM DC bias is adjusted to below the quadrature point 
and the amplitude of RF signal is less than V n 26 . In this way, the pulse duty cycle can be 
less than 50%. In order to induce large frequency chirp to the pulse, an RF amplifier is 
employed to amplify the RF signal driven on PM. Then a segment of SMF is used to 
compensate the chirp induced by the PM. It is known that as the SMF length 
increases, the FWHM is decreased to the minimum, and then is increased 
subsequently. To get a shortest pulse, we need to optimize the SMF length so that the 
pulse chirp can be totally compensated by the SMF dispersion. In the experiment, by 
changing the input RF power along with the length of the SMF, pulses with different 
FWHMs can be obtained. 

Devices fabrication. We employ an on-chip MDR to observe first-order OD and 
OHT. First we design and fabricate the MDR on an SOI wafer. The thickness of the 
top silicon and the buried oxide layer of the SOI wafer are 340 nm and 2 urn, 
respectively. Then the device layout is transferred to ZEP520A photoresist by E-beam 
lithography (Vistec EBPG5000 + ES). Then the upper silicon layer is etched 
downward for 240 nm to form silicon ridge waveguide and input/output grating 
couplers, through inductively coupled plasma (ICP) etching (Oxford Instruments 
Plasmalab System 100). We use vertical grating coupling method to couple the fiber 
and the silicon MDR, and the grating coupler we designed has a period of 630 nm 
with duty cycle of 56%. The coupling loss of the grating coupler is measured to be 5 dB 
for a single side 27,28 . 
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